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Three-Dimensional Magnetohydrodynamical Accretion Flows into 

Black Holes 
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Abstract 

Outflows and convective motions in accretion flows have been intensively discussed recently in the con- 
text of advection-dominated accretion flow (ADAF) based on two-dimensional (2D) and three-dimensional 
(3D) hydrodynamical simulations. We, however, point that without proper treatments of the disk magnetic 
fields, a major source of viscosity, one can never derive general, firm conclusions concerning the occurrence 
of outflows and convection. We analyzed the 3D MHD numerical simulation data of magnetized accretion 
flows initially threaded by weak toroidal magnetic fields, finding large-scale convective motions dominat- 
ing near the black hole. In contrast, outflows occur only temporarily and are not very significant in our 
simulations. If there grow strong vertical fields somehow, however, formation of bi-polar jets is inevitable. 
It is claimed that radiation could be dominant at the outermost zones of the convective disks because of 
outward energy flow by convection, however, this is no longer the case in convective MHD flows, since 
accretion energy can be released via magnetic reconnection in the inner parts. Such reconncction leads to 
sporadic flare events, thus producing substantial variability in out-going radiation. 
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1. Introduction 

Rapid progress has been made recently in the research 
fleld of black- hole accretion flow and our knowledge about 
flow properties has been remarkably enriched. In addi- 
tion to the standard disk (Shakura, Sunyaev 1973), the 
existence of a distinct type of flow, advection-dominated 
accretion flow (ADAF) has been widely recognized. Al- 
though ADAF was proposed long time ago by Ichimaru 
(1977), its signiflcance is only recently fully appreciated 
(Narayan and Yi 1994, 1995; Abramowicz et al. 1995; 
for a review, see Kato et al. 1998). 

Narayan and Yi (1994) made two important pieces of 
predictions: possible occurrence of convection (in the ra- 
dial direction) and outflows. Since radiative cooling is 
inefficient, entropy of accreting gas should monotonically 
increase towards a black hole (i.e., in the direction of 
gravity), a condition for a convective instability. Convec- 
tion was later conflrmed numerically (e.g. Igumenshchev 
et al. 1996; Stone et al. 1999). The occurrence of out- 
flow is, on the other hand, still open to question; the self- 
similar model points an advection-dominated flow having 
positive Bernoulli parameter. Be, meaning that matter 
is gravitationally unbound and can form outflows (see 



also Blandford, Begelman 1999). However, matter with 
positive Be may not necessarily produce powerful un- 
bound outflows (Abramowicz et al. 2000). Certainly, 
the traditional, radially one-dimensional treatments are 
inappropriate to discuss this issue. 

Through extensive 2D hydrodynamical simulations, 
Igumenshchev and Abramowicz (2000, hereafter lAOO) 
have found different modes of accretion flow realizing in 
advection-dominated regimes, depending on the magni- 
tude of the viscosity parameter, a. They found outflows 
for a moderately large a ('^ 1) and convection and/or 
large-scale circulations for small q;<0.1 (see also Igu- 
menshchev et al. 2000 for 3D calculations). 

However, the magnitude of viscosity is not a free pa- 
rameter in real situations but should be determined in a 
self-consistent way. It is now widely believed that mag- 
netic flelds provide a major source of disk viscosity. If 
so, without proper treatment of disk magnetic fields one 
can never derive solid conclusions concerning flow pat- 
terns. Dynamics of magnetically dominated accretion is 
in question. This prompted us to examine the three- 
dimensional data of the global, MHD disk calculations 
flrst made by Machida et al. (2000), aiming to establish- 
ing a general, new view of accretion flow in advection- 
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Fig. 1.. Density isocontours overlaid with a momentum vector multiplied by r (i.e. prv) on the {'co,z) plane averaged over azirnuthal 
angle ((f): (a) snapshot at t = 10.5, (b) that at t = 19.7, and (c) the time averages from t = 10.5 to t = 20.8, respectively. The 
spacing of density isocontour is A logp = 0.1. The density contours where p < 0.1 are not shown. 



dominated regimes. The results of MHD flow structures 
are presented in section 2. The final section is devoted 
to discussion. 



2. Three-Dimensional Flow Structure 

Machida et al. (2000) calculated how magnetic field 
evolves in a rotating MHD torus with its center at uj — rij 
initially threaded by toroidal {B^) fields under the New- 
tonian potential (Okada et al. f 989) with cylindrical co- 
ordinate (-07,0,2). They solved ideal MHD equations, 
but entropy can increase by shock heating. Since no cool- 
ing is taken into account in the computations, the simu- 
lated disk is of an ADAF type. They assumed the adia- 
batic index, 7 — 5/3, and no heat conduction. The num- 
bers of mesh points are (n^ xn^ xuz) — (200 x64x 240). 
The outer boundaries at zu — 6.4 vq and at z = 11.8 tq 
are free boundaries at which waves can transmit. A peri- 
odic boundary condition is imposed for (/)-direction. The 
grid size is Azu — Az — 0.01 ro for < tu < 1.2 tq and 
< z < To and otherwise increases with zu and z. We im- 
posed absorbing boundary condition at r = 0.1 ro where 
r = (n7^-|-z^)^'^. To initiate non-axisymmetric evolution, 
small-amplitude, random perturbations are imposed at 
t = for azimuthal velocity. Although the initial torus 
has a fiat angular-momentum distribution, the angular- 
momentum profile soon approaches that of the Keplerian 
(~ vj^^"^) after several rotations. Since it is difficult to 
add gas threaded by magnetic fields, we assume no mass 
input. However, the fiow can be regarded as being quasi- 
stationary inside ro, since the total calculation time is 
t ~ 22 (in a unit of the rotation period at ro) which 
by far exceeds the mass-fiow timescale (^ r/|wr|) in the 
inner zone. 

As time goes on, each component of magnetic fields is 
rapidly amplified via a number of MHD instabilities to- 



gether with the differential rotation. The maximum field 
strength is determined either by field dissipation by re- 
connection or field escape from accretion fiow via Parker 
instability. As a result, the plasma /3 (= Pgas/Pmag, the 
ratio of gas pressure to magnetic pressure) finally reaches 
~ 2 — 10 on average, irrespective of the initial /3. The 
corresponding a value is a ^ 0.01 — 0.1. Even lower-/3 
(< 1) values are observed in the regions with filamentary 
shapes. Although the spatial distribution of the fields 
is quite inhomogeneous (Kawaguchi et al. 2000), matter 
distribution is somewhat smoother. In the following, we 
analyze the case with initial /3o = 100 at {zu, z) = (ro, 0). 

Figure 1 displays the snapshots (left and middle) and 
the time averages (right) of density contours and momen- 
tum vectors on the {za, z) plane. Note that each quantity 
is averaged over azimuthal angles. Large-scale convective 
motions, which look very similar to those of lAOO (their 
figure 15) are evident in these panels; that is, the accre- 
tion is suppressed in the equatorial plane and the mass 
inflows concentrate mainly along the upper and lower sur- 
faces of the torus-like accretion disk. The polar inflows 
are highly supersonic, although accompanying mass flux 
is small because of low density. Note that we assume 
equatorial symmetry, thus large-scale motion penetrat- 
ing the equatorial plane is artificially suppressed in our 
simulations, while it is allowed in lAOO. 

In figure 2, we plot various physical quantities averaged 
over a spherical shell with a constant r: 
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The results are summarized in Table 1, in which we also 
compare radial dependences of various quantities for var- 
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Fig. 2.. Various physical quantities averaged over angles as functions of radius (r): matter density (upper-left), temperature (up- 
per-middle), magnetic energy (upper-right), radial velocity Vr (lower-left), azimuthal velocity (lower-middle), and viscosity parameter, 
a = —{BrB^/iTr)/ {pg^s) (lower-right), respectively, at different epochs, (a) at t = 10.5 and (b) at t = 19, 7, respectively. The dashed 



curves represent power-law relations, such as p (x r "■^, T oc r 



B-' (xr 



^■'\ and vj, (X r ^■''. Long dashed curves : 



the lower-left panels indicate portions of negative Vr (inflow) . The unit of velocity is the Keplcrian rotation speed at t-q . 



ions types of accretion. Note that outer parts of the sim- 
ulated disks are formed by out-going motions of the disk 
material as a result of re-distribution of the initial disk 
angular momenta as is indicated by positive radial ve- 
locities. In the MHD flow, we concentrate on the flow 
structure at r < tq, where inflow motions dominate (in- 
dicated by long-dashed lines in the Vr distribution). 

In general, our numerical results remarkably agree well 
with those of the convective disk at least at t — 10.5, 
confirming that convection dominates the flow dynam- 
ics. It is interesting to note that density profiles differ 
at t — 10.5 (named as early stage) and at i = 19.7 
(late stage). We notice that entropy [oc \n{pgas/p^) 
with 7 = 5/3] increases with decreasing r in the early 
stage (namely, the flow is certainly convectively unsta- 
ble), whereas it is nearly fiat (i.e., the flow is isentropic) 
in the later stage. This indicates that convection effi- 



ciently works in the early stage to smooth out entropy 
profile, leading to a flat entropy profile in the late stage. 
The plasma (3 is roughly constant radially, whereas elec- 
tron drift velocity (ex j/p oc |V x Bj/p), a key factor to 
turn on anomalous resistivity, increases with a decreas- 
ing r as oc r~^ . Magnetic reconncction, hence, preferably 
occurs at small radii. We also notice that a has a weak 
radial dependence; roughly, a oc r'^-^. 

3. Discussion 

Magnetic fields play many crucial roles in activating 
accretion disks: source of viscosity (Shakura, Sunyaev 
1973), producing rapid fiuctuations (e.g. Kawaguchi et 
al. 2000), jet formation (e.g. Matsumoto 1999), and so 
on. Nevertheless, our knowledge about the disk mag- 
netic fields still remains poor mainly because of their 
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highly nonlinear and global nature, requiring sophisti- 
cated multi-dimensional simulations which are only re- 
cently available. A number of questions have been ad- 
dressed: is the field strength strong enough to account 
for observed rapid accretion? Does convection really oc- 
cur in a magnetized disk? Or is it more likely for an 
ADAF to produce outflows/jets? 

Although we have not run many models, we obtain a 
sort of answers for some cases. The field strength can in- 
deed grow to be strong enough explain observed accretion 
phenomena. Convective motions are evident in the simu- 
lations, indicating that convection may be rather general 
phenomena in accretion fiows. The estimated viscosity 
parameter is, a ~ 0.01 — 0.1, for which we indeed expect 
the onset of convection from 2D/3D hydrodynamical sim- 
ulations. More precisely, we find a oc r^-^ . However, we 
should keep in mind that the presence of patchy low-/3 
regions implies that simple a-type prescription may not 
perfectly work. 

As Machida et al. (2000) already noted, outflow is 
temporarily observed, although it is not very significant; 
the mass outfiow rate is less than the mass accretion 
rate by some factor. Although we began calculations 
with weak toroidal fields, if we start with strong vertical 
fields, magnetically-driven bi-polar jets are evident (e.g. 
Matsumoto 1999). The flow pattern of MHD disks may 
depend on the initial field configurations. On the other 
hand, injection of a magnetized, differentially rotating 
buoyant element naturally can drive outflows (Turner et 

al. 1992] Wn nnod fiirthpr A/THD t;tnHins hpfnrn fir 
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are brought outward by convective gas motions. That is, 
the computed spectra on the basis of the ADAF model 
will not need significant changes even in the presence of 
large-scale circulations. 

We have not yet included conductive energy trans- 
port, although two-dimensional hydrodynamical simula- 
tions found significant changes in flow structure caused 
by conduction; any conductive models did not show bipo- 
lar outflows nor convection, since the thermal conduction 
mainly acts as a cooling agent (lAOO). It is tempting to 
check if this is the case in MHD disks, as well. 
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It is suggested (Igumenshchev 2000) that since energy 
can be carried outward by convective motions, radiation 
could be dominant at the outermost zones. We, however, 
point out that magnetic reconnection is likely to occur in 
the inner zones and to release substantial fraction of the 
total energy of accretion flows in the inner zone. This 
is a very likely possibility, if indeed magnetic flares are 
responsible for variability (e.g. Kawaguchi et al. 2000). 
In fact, we find large field energy in the inner zone, which 
can be released in a sufficiently short time before fields 



Table 1. Various types of accretion flow. 



accretion mode p{r) T{r) Vr{r) B (r) 

ADAF 
outflow 
convection 
MHD flow 
(at t = 10.5) 
MHD flow 
(at t = 19.7) 
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